INSTITUTE  OF  MATHEMATICAL  SCIENCES 

LIBRARY 

?5  Wavcfly  Place,  New  York  3,  N.  ^ 

:L  JULY  4,  1943 


A  SCHEME  FOR  THE  COMPUTATION  OP 

SHOCK  7;a^/es  in  gases  ahd  fluids 

BY   HER?.'IANN  VJEYL 


AM?  MEMORANDUM  33.7 
AMG-KYU  I.-IEMORAKDUM  18 


APPLIED  MATHEI'ilATICS   .^AKEL 
NATIOML  DEFEUSE   RESEARCH    COMMITTEE 


THIS  DOCUMENT  CONTAINS  INFORMATION  AFFECTING  THE  COPY    NO    5^^ 

NATIONAL  DEFENSE  OF  THE  UNITED  SfAIi-S  WITHIN  TM?  MEAN1N8  — ^ 

OF  THE  tSHOMAOE  ACT,  U  B.C.  5"!  M  A""  ^2.  I'  T:A"'3- 
MISSION  OR  Tttt  BFvilA''- N  i.f  I  V- C  '  T*  .-)  l.-i  'V..\  'U-R 
TO  AH    UhAUIMOKIZtO    htliSuN    IS    miyul..lT.i>   UY   i.l>W. 


A   SCKE^.'ui   FOR  T'lE    COI'PUTATIOII   OF   SHOCK  WAVES 
IN   GASES  AITD  FLUIDS 

By 

Hermann  Vieyl 


Temperatures  in  shock  waves  mount  so  high  that 
the  assumption  of  constant  heat,  or,  v/hat  amounts  to  the 
same,  of  the  linearity  of  the  function  f  v/hich  expresses 
the  energy  U  per  mol  in  terms  of  temperature,  is  only  a 
rough  approximation.   Instead  we  are  going  to  treat  f   as 
an  empirically  given  function.   In  order  to  obtain  a  clear- 
cut  algebraic  algorithm  one  is  then  forced  to  express 
everything  in  terms  of  the  temperature.   It  turns  out  that 
ideal  gases  end  fluids  v;hich  obey  Tammann '  s  equation  of 
state  can  be  tackled  by  exactly  the  same  scheme. 

A.  Onc-dimensional_  Shock 
^■1.  Ideal  ^as.  Let  v,  p,  T  denote  velocity 

LI , il. ^  -x 

pressure,  absolute  temperature,  and  V,  U,  i  =  U  +  pV,  ->-^ 

volume,  internal  energy,  heat  content,  entropy  per  mol. 
/^   =  V~   is  the  density.   We  choose  a  coordinate   x   so 
that  the  shock  front  rests  at  x  =  0.,   The  problem  is  then 
stationary.   The  indices  1  and  2  refer  to  the  states  of  the 
gas  before  and  after  the  shock,   x  <  0  and  x  >  0  respec- 
tively.  The  equations  relating  the  states  1  and  2   are 
as  follows 

a)  f^^   ^/gVg  (=  M)  , 

b)  Mv^  "^  Pi  "^  ^''*^2  "*"  ^2   * 
1.2  -  .   ^1,2 


°)      H  "^  2^  =  ^2  ^  t  ^2 
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The  law  of  increasing  entropy  requires 

d-)     lUn^  -  Hi)  >  0   • 
We  determine  the  direction  of  the  positive   x-axis  by 
M  >  0.   Then  the  shock  front  has  the  nog^.tive  velocity 
-  V,   with  respect  to  the  gas   in  the  state  (1),  i.e. 
the  front  moves  into  (l),  or  (1)  is  Indeed  the  initial 
state  in  which  the  gas  finds  itself  before  it  is  reached 
by  the  shock.   Inequality  d---)  now  reads 

^)       ^Y']2  -   ^1>°- 

For  an  ideal  gas  we  have  the  relation  p  =  RTy<^  , 
and  U  is  a  monotone  increasing  function  R.f(T)  of  T 
alone.   Therefore 

i  =  R.j6{T)    ,      j5{T)    =  f(T)  +  T   . 

In  b)   we  replace   p  by  RT p   and  then  use   a)   to 

M 
substitue  •-   for  p  ,   Supposing  the   tamneratures   T, 

and   Tp   to  be  given,  wc  obtain  the  following  two  equations 

for   V,   and  Vp . 

RT,         RT^ 

(2)  R(^(T2)  -  j6{t^))    =   i(v^  -  v|)   . 

It  is  clear  that  they  detemine   v, ,  v   except  for  a 

common  sign.   Honce  we  shall  obtain  unique  values  for  the 

/  2  2  r— 

ratio     v,/Vp      and    the    squfires      v'    ,    Vg    .      By  writing     v^R 

instead   of      v     we   effect   the   normalization     R  =   1. 

Multiply    (1), 

_   ^2        '^l 
-1    -   -2   -  V;   -    v7 
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3. 
by  v^  +  Vp   and  use  (2)  ,' 

Introducing  the  function 

g(T)  =  2j25(T)  -  T  =  2f(T)  +  T 
we  get  for  the  quantities 

"^1    ^1  v^  '  ^2        ^2   Vg 

the  equations 

(3)  ;  1^2  ■  ^1  "  ^^"^2^  "  S(T^),      D^^Dg  =  T^Tg   . 

To  fix  the  ideas  suppose  '^i    "^  "^o    '         ^^^  "^  nothing  hut 
the  simple   q^;  drctic  problem  to  trf-,nsform  a  rectangle 
of  the  sides   T, ,  Tp   into  one  of  equal  area  whose  sides 
have  the  increased  difference 

giTg)  -  g(T^)  -  (Tg  -  T^)  I-  2{f{T^)    -  f(T^))  >  Tg  -  T^  . 

Hence   D  <  T  <  Tg  <  Dg  .   (Were   Tg  <  T^  wo  should  have 
Dp  <  Tp  <  T^  <  D,   ,)   With  the  "displaced  temperatures" 
D. ,    Dp  wo  find 

!  ,. -1 T  ^2  '_  ?i"W  "^ 

(4)  I  1     j2     1 

I  Ps  "  ^i  '  Pl  "  ^1     . 

and  after  some  elementary  trans  formations  the  equation  (2) 

,  T,.  V,    2        D„ 

and  ^i    ^\      =  _^ 

^1^2^      ~\ 

•  2    -      Sli  ^   1  f  2    '-_^1  ' 

^1   "   T;    -    D     '   "^1         2    flTo)    -    f(T7T      ' 
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r= 

^2 

?2 

^2 

• 
9 

• 

r 


2  2 

T      -   D 


In   this  way   all   relevant  qr;antitles   are   e:z:pressed    in   terms 


of  T. 


thus 


ind   T 


2  • 


The  entropy  h  is  p:iven  by 

n/R  =  1   ^^  qT  -  log  O  , 


1(^2-  1l^  =  Ji^-^dT-  lo^ 


where 


R 


V\ 


T  ) 


T  T 


■!^  I 


D, 


ir^dT  +  loe:nf^  . 
i  ^1 


This  function   h(T,  ,Tg)   is  clearly  skevz-symmetrlc.   We 
make  the  hypothesis  which  is  fulfilled  for  most  if  not  all 


gases  in  nature  : 


(K^): 


I  n  (T^jTg)   has  the  sign  of  T^  -  T^   . 
For  a  small  difference   of  T,  ,  T   an  easy  calculation 


yields  , 


'I  ^^l'^2^  "  F'^^'^^-Vt")   ^'^^^  T^  =  T  -  I  /\  T,  Tg  =  T  +  |A  T, 


Wf  =i[f'(f'  +  l)(2f'  +  1)  -  f"T]  . 


where 
(5) 

Hence  our  assumption  of  the  monotony  of   f  and   (H, )   imply 
Df  =  f '  >  0,        Wf  >  0. 
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Question:   Do  these  differential  relations  vice  versa 
imply  (H  )  ?  i7e  shall  see  that  the  expression  Wf  occurs 
again  and  again  in  our  investigation. 

Notice  that  A*]  is  of  the  order  (/3.  T)"^  for  a 
small  difference   ,.A,  T  =  T^  -  T,  .   On  the  basis   of  (H^) 
the  inequality  d)  is  equivalent  to   T,,  >  T,  :  the  shock 
raises  tempera turo  and  density  (condensation  shock). 
$2.  Fluid.   In  Tammann's  equation  of  state 

(p  +  K)(V  -  b)  =  CT        (b,K, C   constants) 
wo  introduce 

p'  =  p  +  K,      V  =  V  -  b. 
The  fundamental  thermodynamical  equation  for  virtual 
variations 

T^v   -   ^U  -r  p  J  V 

gives 

<^  \  — t"     ^p   "  ^' — f T  ^nr*^  ' 

Consequently     ^ — "^-^ — -     is   a   total  differential,    or 
U  -  KV    is   a   function      C.f(T)    of     T      alone: 

U  =  IW    +    C.f(T), 

1   =    (pV    +   bp)    +  KV    +    C.f(T) 

=  bp   +  p'V»    +    Cf(T)    =  bp   +   C.jZ5(T)    . 

Mv  +  p  has  the  same  value  for  the  states  1  and  2,  and  so 

1   p 
has   1  "•■  "I  V  .   By  subtracting  b   times  the  first  expression 

fi-om  the  second,  one  gets 

C,J2)(T)  +  "I  v^  -  Mbv  , 
or  aftor  replacing  v  by  l^'     on  the  basis  of   a), 
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one  then  finds  the  following  two  quantities 


cross  the  shock  front  unchanged.   Hence  the  sane  is  true  for 


M^'V  +  p'  =  M^'V  +  -^W  and 


Setting  MV  =  v'   we  get  the  same  equation  (1)  and  (2) 

as  before : 

CT         CT 


After   V ' ,  v'   are  determined  from  them,  the  relati 


on 


T'» 


T/!  =  V  '  /V  = 

^'-    V  1    V|  -  b 

gives  the  value  of  M   if  the  3:wCific  volume  V,   of  the 
fluid  in  its  initial  state  1  is  known  besides  the  tem.perature 
T^.   Then  Vg,  v^,  v^   follow  from 

V   -  b   vj;  V 

v^  -  v^  =  V'  -  v^   - 
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B.  Reflection  of  _a  Shock  '/Vavo 
1.  The  formulas.   In  this  tv;o-d.imensional 

problem  we  have  two  shock 
fronts  1  and  2  at  rest,  which 
X  meet  on  the  reflecting  wall. 

-^^'-.  -i^"^  t".      The  components   (x,y)   and 

■"■"---^     -'   ''-'''  '  t  L 

«^^>-^-. /'^ ''•'■   ^    I.  . ..  (x,  ,y.  )   of  any  vector  in  the 

(' 

coordinate  systems  L  and  L, 
(see  figure)  are  connected  by 

x^  =  a^x  -  b-j^y,     x  =  a^x^  +  b^^y-j^  ,  /a^  =  cos^x^"'  , 

!  *         ) 

y^  =  b^x  +  ap-      y  =  -'b^x^  +  a^y^   \b^  =  sin  oc^  / 

Similar  relations  connect  L  and  Lr-,    .   Let  (u-,,0),  (u,v), 
(Up,0)   be  the  com-poncnts  v:ith  respect  to  L  of  the 
velocities  in  the  three  regions  marked  by  their  temperatures 
T  <  T  <  T^  .   The  components  of  the  first  tv>'0  velocities 
with  respect  to   L,   are 

The  tangential  component  has  the  same  value  on  both  sides 
of  1,  hence 

1  ^1 

^1  =  I^'^V  -  V^  '     Vi  =  TT^^ V  -  V^    • 

Thus  the  quotient  of  the  normal  velocity  components  in  (T^) 

and   (T)   has  the  value 

b^   a^u  -  b^v 
"1   a^  b^u  ^  a^v 

Again  we  are  going  to  express  everything,  including  the 
angles  oc    c<-„   of  the  incident  and  the  roflcctec  waves, 
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in  terms  of  the  tenperatures,  thi;s  reducing  the  problem  to 
simple  universal  algebraic  equations  once  the  energy  function 
U  =  R.f(T)   is  known.   (One  mi^lit  be  inclined  to  consider 
T^,  T,  C/ ^   as  the  "natural"  data,  but  then  the  problem 
is  not  algebraic.)   The   v   in  a)  and  b)  is  now  the 
normal  velocity  component  while  v   in   c)  is  the  square 
of  the  total  velocity.   But  because  of  the  continuity  of 
the  tangential  component,  even  in   c)  the   v  may  be 
interpreted  as  th.e  normal  component.   Hence  the  same 
relations  (3)  and  (4)  as  before  result  for  both  fronts  1 
and  2: 


D'  -  D^  =  g(T)  -  g(T^)  , 
D"  -  Dg  =  g(T)  -  gCTg)  , 


o 


"^1='/- 


2_  =  iil 


"1 
D" 


> 

D'D^ 

= 

TT^ 

• 

i 

D"D, 

= 

TTg 

• 

(.- 

=  )^^ 

/' 

= 

^1 
^1 

i 

ll 

^1 

m 

P 

T 

> 

c-i' 

=•? 

= 

> 

£2 
P 

- 

T 

. 

(6) 


P2  "  ^2  ' 


For  the  velocity   (u,v)  we  have  the  two  equations 

a,u  -  b-.v    a 

b"u  +  a^v    bT  •  ^1 


(b^u  +  a^vf   -  [~^(a^u  -  h^v)f   =  2{jiiT^)     -  i?i(T)). 


The  first  equation  is  equivalent  to  the  relations 


a^u  -  b^v  =  -'^^  .  a^k^ 


b.u 

X 


a,v  =  X 


.  b. 
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involvinf3   an   unknov/n  factor    ^- ->  •      They   CF.n  bo    solved  v/lth 
respect   to     u     and     v     as   follo".vs: 

(8)  u  =  a^    .  A^a^V:^   +  b^    .  A  ^b^   =    '>^-^i^l   +  a^k^)      , 

(3')    V  =  -b^  .  ^;;^a^k^  +  a^  .  X^b^  =  X^a^b^d  -  k^)  . 
The  left  member  of  (7)  is 

•^[a^(b^u  +  a-|V)  +  b^(a -u  -  b^v)  ]  = 
^1 

^(a^  .  X^b^  +  b^  .  X^a^k^^)  =  ^  .  '^-^a^b^d  +  k^)  , 
^1  ^1 

ai  d  if  one  dotepmines  the  factor  X^  by  (8'),  the  equation 

(7)  turns  into 

(9)  ^  .  -. ^  =  2[/(T^)  -  ^(T)]   , 

a  J    ^   "    '-1        ^ 

while    (S)    and    (3')    give 

a  b,  { 1   -   k,  ) 

(10)  y  _      X   X  1 


Equations    similar   to    (9)    and    (10)    prevail   for 
the   second  front  2.      Comparison   leads   to   the   follov;lng 
two   relations   for   the    tv/o  unknown   angles: 

^A^^    -   ^1^      =   ^2"?M   '    ^2^       (^  v^ 
^1    "   4^  4   "    4^2 

Vi  =  V2    ^=^^^' 

Here  v;e  have    set 

A^   -  2(^(T^)    -   ^(T))    =    (g(T^)    -   g(T)    )    -;     (T^    ••    T), 
1   -   k. 


^1   ~  -^11   +   k,       ' 
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while  A^,  B^  have  the  corresponding  meaning.  Use  the 

tangents  _1   _2   as  the  unknowns   t, ,  t^  : 
^1  '  ^2 

t^(l  -  k-^)  _  tgd  -  kg) 

^1  "  ^1        ^2  "  4 

B         B 

^  -    ^     (=  s)  . 


2         2 
1  +  t^    1  "^  tg 

The  last  equation  gives 

2   ^1  2   ®? 

h    s    ^'  ^2   s    ^   ' 

and   then,    after   squaring  both  sides   of   the   first  equation 

one  gets 

B,    -   s  B     -    s 


v/here 


(C^    -    s)^         (Cg   -    s)2 

^1  =  r^  =  -i-TT^  =  ¥h  ^  ^i)  '  S  ^  h^2  ^  \^ 

The  resulting  equation  for   s   , 
s  -  Co      s  -  E., 

is  not  cubic  as  one  might  fear,  but  only  quadratic.   Setting 


s    -   B^        ^ 

we   find 

z   -    1      ' 

^-    ^2 

■     s    >    C^ 

z  -(q^z  +  qg) 
1  -(q-i  z  +  q^) 

1.  i^ 


where 


^1  "  ^1  ^  ^      ^2  ~  ^  ^  ^ 

c^  -  Cg-^i  >  Cg  --C,  -% 
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The   equation    (11)    becomes 

.  z    -    (q,  z   +   q^) 

1  o 


-.2 


>  1    -Tq-|_z   +  'qgT/ 


-    1   =   z    -    1 


nd  fter  the  removal  of  the  foreign  factor   z  -  1  : 


(q^^z  +  Qg)   =  z 


or,  finally 
(12) 


q^z^  -  (1  -  2q^q2)z  +  q^  =  0. 


Exhibit   the    inherent   Sj-mmetry  by   setting      z      =   z,    z„   =   l/z : 

J.       ^ 

(12'.)      q^z^  -  (1  -  2q^qg)  +  q^z^  =  0;   z^Zg  =  1   . 

Before  going;  on,  express  the  relevant  quant:' ties 
of  index  1  in  terms  of  the  temperatures  T^ ,  T  and  the 
displaced  temperatures  D^ ,  D   ;  one  readily  finds 

I     A^  =  (L^  -  D')  r    (T^  _  T) 


B, 


(D^  +  D 


')  -  (T^  +  T) 


D'  -  T, 


C,  =  B,  -  T 


Qi  = 


1   Do 


D. 


Equation  (12)  is  simplified  even  more  by  introducing  y z  =  i^    : 


/ 


s  -  Cg 

n"  s  -  c'^  ' 


,  L 


2 


S  -  C, 
s  -  C, 


It  then  reads 
(13)  q 


i^-^i  ""  ^2 ^"^2  =  i>    ^n^ 


1  . 


One  easily  sees  that  the  sign  is  correctly  determined,  namely 
so  that   ^1^1''"  ^2^--  ^^^    *^'''  "^^'^^"^^  """l  ^"'^  ^^'^      '1* 

After   z   or  C-,  has  been  determined  one  finds 


z  -  1 


r.  1 "  (!S^  2 
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'^B^   -   B. 


t,    =   tan  Cx,    =  !-2~ 


12. 

,1/2 


Jh    -   ^2        \ 


\l/2 


^l^l  -  ^2/ 


'1    ""''  ^  1 

For   the   latter  relations  we  might  v;rlte 


,      -tg   =   tanag    -^  ^  ^ 


\  2   2 


1/ 


-    tg/t^         .-^.^,  -    t^/tg   - 


-/ 


Bi      -      B^ 
t,  t^   =      ^ 


a   2        B^6i    -   BgC,^      • 


In  order   to  express     v/u     we   observe    that  according 
to   our   above    calculations 


t^(l   -  k^)    ^  B^   -   s 


and 


Bo  -   B^  p  C^  -    C, 

^1   -   ^ 


i=l   -    z    (=   1   -C'),       c^ 


1    + 


therefore 


t-d    -   k^)    _   Eg    -   B^ 


md   as 


^I'h 


^2'   \  _ 

c     -  r 
^2       ^1 


Cg  -   c^    •    1  -    ,     > 


1    -    q;^    -    q. 


pnd    (1   - 


r  \,-    - 


I  '' 


=   t^    +   t„,    the   e:?pression 


t^(l  -  k^) 
__ 

t   +  k 
^1    ^^1 


for  v/u  may  be  written  in  the  more  symmetric  form 


V 

u 


^  -  ^1  -  ^2 

^1  ^  ^2 


g  2.  Discussion  of  signs.   The  qualitative 
discussion  of  the  results  is  more  difficult  than  the 
straightforward  algebraic  part.   B  ,  B^,  q, ,  qp   are 
positive,   Cg  >  0,  C^  <  0.   A  first  condition  for  the 
existence  of  a  solution  corresponding  to  the  temperatures 
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T, ,  T,  Tg   is  the  reality  of  the  roots  of  the  equation  (12) 
or  (13),  namely 

(15)  <4  =  1  -  4q^q2  ^  °* 

But  this  is  not  quite  enough;  for  at  least  one  of  the  roots 

z  =  z',  z"  of  (12)  the  quotient   "2  "   1   must  be  positive, 

B^z  -  Bg 

Exactly  one  of  the  quantities  B,z'  -  Bg,  B,z"  -  Bp  will 
have  t&e  sign  of  Bp  <^  B^  ,  if  their  signs  arc  opposite,  or  if 

(16)  q^(B^z'  -  Bg^B^z"  -  B^)  <  0. 
The  left  side  equals 

q^[B^z'z"  -  B^B2(z'  +  z")  +  3^]=  B^q^  -  V^d  -  2q^q^)    +  B^q^  , 
and  hence  the  relation  (16)  reduces  to 

(17)  I    J  =5  (B^qg  +  Bgq^)^  -  B^B^  <  0   . 


If   z',  z"  are  conjugate  complex,  (15)  is  necessarily 
positive.   Consequently  (17)  implies  (15);  the  relation 

(18)  (B^qg  +  Egq^)"  -  (B^qg  -  \^^)^   =   4qiq2B;,_B2 

confirms  this  in  a  more  explicit  v/ay.   Thus  the  case  of 
one  solution  is  characterized  by  J  <  0. 
If,  hoTvover 

(19)  J  >  0,  Q  >  0 

we  may  have  no  or  tv;o  solutions.   The  left  side  of  (16)  is 

nov;  positive,  therefore 

either     B,  z '    <  B^   ,    B,  z"    <  B^  and   thus   E^z'z"   <  B^ 

1  i;  1  <s  1  d 


^7  «  7"    >   R^ 


or  B,  z'   >  B^    ,    B   z"   >  B^  and   thus   B^z'z"  >  B 


The  alternative  takes  on  the  form 

CONFIDENTIAL 


14. 

Result:   In  case  (19)  vie   have  two  solutions  or  no  solutions 
according  to  vih ether 

^1^2  -  ^2^1    ^"^   ^1  -  ^2 
have  equal  or  opposite  signs.   We  add  the  observation  thit, 

due  to  the  relation  (17)  or 

(B^qg  -  Bgq^)^  =  J  +  QB^Bg 

B  q   -  B  q   can  not  change  sign  in  an  interval  whore  (13) 
hold?.  Moreov:.r  (14)  shovi's  that  B^  >  B,  implies   q,  +  qp  <  1 

and  a  fortiori  Q,  >  0. 

^1   ^2 
The  difference — ^  ,  the  sign  of  v/hich  plays 

a  role  in  this  discussion,  is,  by  definition 

T  -  D^     D^  -  T 
(Dg  -  l;^)[2  -  ^T^r^^   -   Tg  -  D"^   • 

^""^  T  _  ^1     D-   - 

Di"  ~  T^  '    Tg    D 

hence  yjo    can  write  more    simply 

and    the   alternative        1  2      ^     ^     amounts    to 


{ 

We  proceed  to  evclurto  this  general  analysis. 
Head-on  collision  is  characterized  by  B^  =  B^    ,    or 
q-  +  q  =  1,  whcih  gives   J  =  0.   Grazing  v/ave  ( t,  =-:x~-) 
is  indicated  by  E,z  -  Eg  =  0,  or   z  =  Bg/B^  must  satisfy 
the  equation  (12),  which  again  gives   J  =  0.   Hence  J  should 

CONFIDENTIAL 


15. 


split  into  two  factors,  one  being  B,  -  Bp.   Indeed, using 
the  expression 


one  easily  gets  p      p 

(B,  -B.)(B^C^-BgC^) 

V  ——————————————  ^ 

Hence  the  sign  of  J  depends  on  those  of  the  differences 

4        ^? 

For   a    .?iven      T     we   plot     Bp  =  B{Tp)    and     ~     =1    (T    ) 

^  ^  -2       '^ 

as  functions  of   T^.   For  small   Tp  -  T  =  „\  T  one  finds 

2 

B(T2).^-f'(f'  +  1)  .^^AP-      ,      V   (Tg)  --.i^  T  -f  Wf(T)  -^^ 

Notice  the  differential  expression  Wf   in  the  formui^a 

for  r'(T2)  -r(T).   Th     kes  it  plausible  that  F(Tg) 

is  a  monotone  increasing  function  of   Tp,  or  at  least  that 

(Kg):  rCTg)  -  r  (T)    has  the  sign  of  Tg  -  T. 


The  behavior  at  the  ends   Tp  -  0  and  Tp  — *- oo   supports 
our  belief  in  the  monotony  of  Y^'lTp);  indeed  for   small  Tp 

4    ^1 

and  I  --»  c<'  with  T^ -^>  cO  .   (H^)  adnitted,  r^  -  -^     will 

alv/ays  be  positive  for  T,  <  T  <  To  »  and  a  grazing,  or 
nearly  .n:ra2lng  ways  proves  impossible .   The  alternative 
J  ^  0  reduces  to  the  sim.pler  Bp  ^  B  , 

B(Tp)  behaves  in  the  neighborhood  of  Tp  =  T 
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according  to  the  formula 

3(T2)  /-^^  f '(f>  +  1)  .-^-^   . 
One  will  therefore  expect  that  the  graph  of  B(T^)  has  the 
shape  of  a  trough.   B(0)  =  2f(T)   and  B(Tg)   tends  to 
infinity  with  Tq^"^  ^"^  >    under  the  hypothesis  that 

f(T)-^00  for  T^(r<i  . 
/    Given  T,,  T,  the 
rnV^/^     /'      horizontal  through  the 
^  /  point  (T   R(T  ))  on 

^^,.--''1"         y        ^  ^^      the  graph  of  B(T2) 
B^r^-i'^" — i ■— "-.  determines  a  value 


T' 


Tg  >  T  for  which  B(T2) 


=  B(Tt)  =  B^ .   For  all   T^  >  T°   the  condition  B^  >  B, 
is  fulfilled  which  indicates  that  the  problem  hcs  one 
solutj  on. 

In  the  interval  T  <  Tg  <  Tg  where   Bg  <  B^  , 
J  >  0,  one  has  to  discuss  the  signs  of 

;d      Dg\, 

Q,  =  1  -  4q^q2      aid  J.  ^   2   -\-~  +  -^ 

Q  >  0  and  L  <  0  being  the  criteria  for  two  solutions. ■  Let 
us  therefore  investigate   1  -  4q,qg   or 

(2/?)  (Dg  -  D^)^  -  4(D'  -  T^)(T2  -  D" ) 

for  T^,  Tg   near   T.   Setting  T^  =  T(l  -  x),  Tg  =  T(l  +  y) 

one  obtains  in  first  approximation 

T^(f'  +  l)^(x  -  y)^  (-  0). 
V'/ere  the  quadratic  form  on  the  left  positive-definite,  one 
would  expect  Q   to  be  pos  itive  throughout.   But  it  is 

COi^IDENTIAL 


17, 


only  ser.l-def inite,  with  the  neighborhood  of  the  diagonal 

2t  2 

X  =  y     as   the   dcaifrer   zone.      If   one    sets     y  -   x  =  ^,-^  ■■;  »x 

p  4 

with  an  arbitrary  coefficient   t  ,  x-h-m  (20)  becomes  T  x 

times  a  quadratic  form  of   t»  the  discriminant  of  which 

once  more  involves  our  familiar  V/f,  namely  equals 

2(1  +  f')1/V  f.   Hence  this  form  is  indefinite  and  Q  becomes 

actually  n-~;gative  in  the  neighborhood  of  the  diagonal 

X  =  y.   The  values  of  Q  for  Tp  =  T  ,  T,  cc-  are  cx^,  1,1, 

Q   I  \  respectively,  and  for 


\ 


given  T- ,  T  we  exT3ect 
1' 

a  graph  of  Q(Tg)  as  shown 
^-^-    in  the  figure.  The  inter- 
val (T^2*  Tjg)  in  which 
"t"    T  t7^--'''t;   XT   "^2   Q,  <  0  is  pert  of  (T,T°). 

L(Tp)  has  a  constant  sign  throughout  the  piece 
(T,  T^g),  ^^-nd  since  ^Z'^.    <   1   and  D^/Tp  =  1   for  Tg  =  T, 
it  is  positive  there  and  hence  no  solution  exists.   For 


-i^ 


T'   <  T^  <  T^  one  will  have  two  solutions,  if  the  constant 

J.  <^        ^        ^ 


sign  of   L   in  this  interval  is  - 


Thus  the  question 


arises  whether   L   changes  sign  before  Tp   reaches   Tp. 
For  a  given  T^  =  T(l  -  x)  near  T   the  expression   ^C^o^ 

actually  changes  sign  at  a   Tp  =  T(l  +  y)   for  vuhich 

2 
y  -  X   is  of  the  order   x  ,  and   Eg  will  be  less  than  B, 

at  this  T   ,  provided 


(H'"): 


W'"f  =  i[f'(f'  +  l)(2f'  -  1)  i-  f"T]  >  0, 


If  we  make  this  further  hypothesis,   L  turns  out  to  bo 
negative  in  the  interval  T^g  <  Tg  <  Tg   for  sufficiently 
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small   T  -  T^   and  hence,  for  reasons  of  continuity,  will 
keep  this  sign  under  all  circumstances:  there  are  then 
ty?o  solTitions  in  that  interval. 

The  fact  that  L  <  0  or  B^q^  <  0   for  Tr.  =  T° 

(B  =  B- )   implies  that  there   q-,  -  Qp  ^  0   or 

(20a)  (Dg  -  D^)(q^  -  qg)  =  (D'  +D" )  -  (T^  +  Tg)  >  0. 

For  the  same  value  of  Tp  vire  have   q   +  qp  =  1,  and  hence 
U^  =  1,   ^2  =  1;      ^n'=   qg/'ll'    ^-"2  =  ^l/^2 
are  the  tv.-o  solutions  of  (13).   In  view  of  the  genorril 
formulas 
(20b)      (tan  C,^)-   =  c..^^_l— ,  (tanC^^)   ==  ^.,  V'^V^ 

the  second  solution  gives  in  the  neighborhood  '^^f  Tg  =  Tg : 


2 


;  q^  -  qg 


Being  real  for  Tp  <  Tg   it  ceases  ,to  exist  beyond   Tg  • 

I'h^  f i.J' s t  solution  must  therefore  be  the  one  v.'hich  continues 

o  / 

beyond  Tp.   For  an  infinitesimal   B  -  B^   set  '-'-,'--  1  +e 

(€  infinitesimal),  consequently 

^g^^l  -  e    and   q^(l  +7.  )  +  qg(l  -  t  )  =  1   or 

B   -  B 
(q^  -  q2)-  =  1  -  ^1  -  ^12  =  ^r^-^      • 

Thus  one  derives  from  (20b)  for  the  first  solution: 
^.    =   -'^-^  =  ct°   at  T^  =  T°  where 

(f.^.^0,-2   _  (^1  "  Bg)  -  (Dg  -  Di)(qi  -  qg) 
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By  (20a)  and  the  expression  of  B,  and   B^   this  formulr. 
reduces  to 


tanCX°  ='V 


;D'  +  D"  -  T,  -  T, 


D^  +  Dg  -  2T 


(Tg  =  T°) 


If  To  tends  to  infinity  then  q,  -^   0  &nd  (-, 
in  such  a  manner  that 

1  :    ^.^  —  l/qi  ,   ^2'~^^1 


>x> 


^/^ 


2      2 
Hence  t,  and  tp  tend  to  zero  and  infinity  at  least  as 

strongly  as  Tg   and  f(T.p)  respectively.   The  limit  is  the 

strange  phenomenon  of  a  head-on  impact  (Cx  -^o)  that  slides 

off  (o<.p  ->  T^)  instead  of  being  throv;n  back  and  thereby 

acquires  infinite  strength  (Po/p  -   Dq/T  — ^C-^^)  .   The  fact 
that  t,  — ^0  with  To— *- •'-^^'  shows  that  OC,  nevoi'  exceeds  a 


certf^Jn  maximum  value  ■••,   < 


T  • 


The  curves  she;;  the 


'^- 


\ 


/ 


i7 


^.. 


or.  yJ 


CK- 


qualitative  behavior  and  are 
'i  in  accordance  with  v.  Neumann's 
■  I  results  for  constant  specif io 
/  '  heat.   In  Nature  the  Mach  phe- 
nomenon of  detachment  bridges 
the  gap  rj.    <  o;-,  <  ^  ^nd 
cuts  off  the  top  part  of  the 
picture  corresponding  to  tem- 
peratures To   above  a  certain 
;  limit. 

i       Our  analysis  reveals  the 
!  importance  of  the  condition 
^♦fwf  >  0.   The  hypothesis  (H") 
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plays  a  more  accessory  role;  as  a  matter  of  fact,  the 

situation  would  be  simpler  if  V/"f  <  0,  because  then  only 

the  branch  T  >  T°   of  solutions  would  exist, 
2     ''' 

C.  Proof  of  the  hypotheses  ( H  )  and  ( Hg ) 

1,  Qualitative  character  of  the  adiabatic  curves. 

We  propose  to  show  that  the  inequality 

Wf  >  0 

is  equivalent  to  the  fact  that  the  adiabatic  curves  in  the 

V,p  -  diagram  are  convex: 

.,2  V 

i^i      >  0. 

'd^'adab. 
Indeed  this  proposition  is  an  immediate  consequence  of  the 

following  simple  formula; 

(21)  1^=2.  ^^^l^XL.  £^     . 
\dV7^^      (f'(T))'^  v-^ 

Proof,  Along  an  adiabatic  curve 

Rf ' (T) .dT  +  p.dV  =0       or 

(22)  dT  _   p 

dV    '  R.f'iTTT 

Because  of  RT  =  pV   the  same  differential  equation  may  be 
written  as 

f'(T)(p.dV  +  V.dp)  +  p.dV  =0    or 

ro-r\  dp  _  f   +  1   __P  _  p    f   +  1  rn    1 

(23)  -^ J, ^.-  R.  -^rr—  T.  -^   . 

Hence  by  logarithmic  derivation  with  respect   to  V, 
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dfp  /  dp  ^   f "    f "    1  •  dT   2 

dv2  '  dV  '  '^'    -^"1  "  ^'''    ^.  ^  ■  V 
-f"'T  +  f «  (f'  ^_J_)_  dT    2 
=    "'^f'Tf  +  ly-T   'dV  '  V  ' 

or,  after*  inserting  (22)  and  subst:tuting  p  V  for  RT  in 

the  denominator, 

:'f"-T  -  f '  (f '  +  1)  _  2  '^  i   =  -  2   __    'Vf 

I     f'2(ft  +  1)       /  V       V"  f.,2^^,  ^    ^^ 

Multiplication  by  (23)  yields  the  desired  rosult  (21). 

2«  ?roof  of  the  hypothosj ^  (H, ) .   With  a  given 
function  f(T)  one  ct.n  form  the  skev/- symmetric  function 
f(T2)  -  f(T^).   It  hr.s  the  sign  of   at   =  Tg  -  T^  for 
an  infinitesimal  c':^!!      if  and  only  if  f '  >  0.   A  class- 
ical theorem  of  Calculus  strt  os  that  this  condition 
f ' (T)  >  0  implies 

f (Tg)  -  f (T^)  >  0   for  Tg  >■  T^ 

even  for  arbitrary  _f_inite  differences  £T. 

Let  f  (T"*  be  any  monotone-increasing  function  of 
the  positive  argument  T.   We  are  going  to  prove  an  analogous 
theorem  concerned  v/ith  a  much  more  conplicated  skev/- 
sjTTBnetric  function   '■"iC^t  ,  Tg)  formed  by  means  of  f  (T)  , 
namely  ^^  ^^ 

r,(T   T)  =   \  i-4r^dT  -  log  -^  =   \-~^clT  -  log 


where  D,  ,  P    f-re  defined  by  the  relations  (5),  g(T)  = 

T  +  2f (Tl .   It  has  been  seen  that  this  T|   for  infinitesimal 
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At  =  Tg  -  Tj   l3  of  tho  order  {  ^,-^;^     ^na 
sign  of  Zl  T  if  and  only  if 

2*V;f  =  f»(f'  +  l)(2f'  +  1)  -  Tf"  >  0. 


«,  f  ^  ,>o  ■ 


¥e   ma 


In  tain  that  \Vf(T)  >  0  for  all  T  implies  ri  (T^,  Tg)  >  0 


for  all  T-j^,  Tg,  Tg  >  T^ .   The  question  asked  at  the  end  of 
k. ,     f,  1,  is  ansv/ered  thereby  in  the  affirmative. 

The  proof  is  suggested  by  the  study  of  a  linear 
shock  v/ave  with  the  temperatures  T,  ,  Tg.   In  purely  mathe- 
matical disftuise  it  runs  as  follows. 

Since  Vjf(T-,)  >  0  implies  'n(T^,  Tg)  >  0  for 
Tg  >  T,  sufficiently;  near  to  T,  ,  "^he  function  -n  Ci'-,  ,  Tg) 
of  Tg  could  never  becvome  nO/^ative  unless  it  passes  through 
zero  for  some  definite  value   Tg  >  T-,   of  the  argument. 
We  will  show  that  this  relation 

(24)  n(T;^,  Tg)   =   0 

leads  to  a  contradiction. 

Choose  tw3  positive  numbers  p^ ,  V,,  determine  R 
by  PnV,  =  RT,   and,  in  agreement  with  the  scheme  (4), 
Pg,  Vg    by 

P__2_D2    Tg        ^2_£_1    ^2 

Pi  ~  T{   .Di  '    ^1  "r":p'2 
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These   relations   inply 


PpV  T 

■p~Yl   =   T^-      '  h^^^^      P2"2   =      ^Tg. 

The  hypothetic  equation  (24)  may  then  be  written 

Ti 
:.(T,,  To)  =   ■   -V"  ^^  ^  log  :^   =   0. 

Form  the  f --.llov.-ing  "entropy"  function  'M  "      of  p,  V: 

-r,'"'  =    j  — -^"--dx  +  log  I-  where  T  =  ipV , 
■T  -^ 

It  has  been  shov/n  in  :^'l,  that  along  the  "adiabatic  lines" 
(Pj  V)  =  const,  one  has 


n 


^P  ^  o      IIP  >  n 

qV 


(24)  can  now  be  interpreted  as  the  statement  that 

the  two  points  (V,,  p, )  and  (Vp,  p^)  lie  on  the  same 
adlab.-.  tic  T'l"  =  0.   Because  of  Tp   >  T,   we  have  Vg   <  V,. 
Along  the  adiabatic  -n  "  =  0  the  following  differential 
equation  obtains 

f'(T)-cT  H-  -^dV  =  0   or    Rf '  (T)  dT  +  p-dV   =  0 

and  thus  we  fnnd  by  integration 

A 


R[f(T2)  -  f(T^)^   =    \  p  dV 


'2 
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T  he  right  side  is  the  shaded  area  F  buloM.'  the  adir.br'. tic 

Y)  =  0.   However,  there  is  a 

■( 

second  more  direct  way  of  deter- 


.-■  X  ^■'X\ 


mining  the  difference 


Y, 


R(f(T,,)  -  f(T^)).   Indeed,  2(f(T2)  -  f(T^))  =  g(Tg)  -  g(T^)  - 


\ 


-/ 


R(f(T2)  -  f(T^))  =  \   P^V^Cl  +  ^  )(1  -  ^)      or 

^'1        1 

RCfCTg)  -  f(T3_))  =  -^(p-L  +  P2)(V-^  -  Vg) 

(Hugoniot  equation).   1'hus  we  have  found  two  values  for 

R(f(T^)  -  f(T-,))?   the  area  of  F  ai  d  the  areaof  the  trapezoid 

obtained  from  P  by  straightening  the  upper  boundary  21.   The 

contradiction  is  obvious;  because  of  the  convexity  of  the 

curve  the  second- area  is  larger  than  the  first. 

3..'?j.Il$S.^-.?S.''9-lii?y^^^l^    ^^^^'^)  '      ^^  ■'^^  write  T.,  ,  T  2 
instead  of  T^  Tp  the  proposition  (H^)  concerns  the  difference 


^l-*-  ^2  -  ^1 


f'(T-^)  +  1 


-  T 


f'{T^) 


n 


B  ut 


T^(D,  -^  Dg  -  T^-  Tg  )  =   T^D^  -   T^D^  -  T^ 

(Dg  -  T^)(T-^  -  I'3_). 


^^1^2 


Thus,    our   tasV  is    to  prove; 
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Dp  -  T,     f  (TJ  +  1 

has  the  sign  of  Tg  -  T^  ,  provided.  Vvf  >  0. 

This  time  we  are  not  dealing  with  a  skevw symmetric 
function  of  T   T..,  ,    and  we  shall  see  that  the  situation 
for  Tp  >  T,   is  different  from  T.  <  Tg  .   Nevertheless  we 
can  give  a  proof  holding  uniformly  for  both  cases. 

Once  for  all  we  suppose  Wf  >  0.   Considering  that 
our  calculations  in  B.,   §  2,  have  established  our  theorem 
for  Tg  sufficiently  near  to  T,  ,  "t  is  suf-'^icient  to  show  that 
the  equation   V'  (  T,  ,  T.-,)  -  0  can   never  hold.   The  proof 
Is  indirect:   we  assum.o  that  v/e  are  in  possession  of  two 
different  values  T, ,  Tg   for  which 

y(T^,  Tg)  =  0   . 

We  choose  and  determine  the  numbers  p, ,  V  ,  pg,  Vg, 
R  as  before.    Then  (25)  gives 

Pg  -  Pj  _  f ;t^)   +  1   p, 

^1  "  ^2      f  (T^)    V^ 
We  now  introduce  the  straight  line  through  the  point  (V^,  p^^ 


defined   by 

P  - 

Pl 

\- 

V 

f'(T^J  +  1    P3 
f'(T^)     V^ 


26  . 

Because  of  the  equation  (cS)  this  is  the  tangent  of  the 
adiabatic  ''')  "  =  0  at  the  point  (V,,  p, )  .    Our  hypothetical 
equation  (25)  states  t?iat  the   point   (Vg,  Pg)  lies  on  this 
tangent. 

Integrating  the  relation 

.  RT-drf"   =   Rf'(T)-dT+  p.dV 
along  the  segraent  12  of  our  tan.-'ent,  we  get 

,  2 


R  \     f '(T) -dT  =  RCfCTg)  -  f(T^)) 


and 

■L  I 


\  p-dV   =   -  (     p  dV   =    -  |(p^  +  Pg)(V^  -  Vg). 
According  to  the  Hugoniot  equation  these  two  expressions 


cancel  one  another;  therefore 

(26)  \T-d-n   =  0. 

I 

But  becur.-->?-of  the  convexity  of  all  adiabatic  lines,  t'] 

varies  monotonely  (namely  decreases)  when  one  moves  along  the 
tangent  starting  from  1,  whether  in  the  one  direction  or  the 
other.   This  circumstance  renders  the  equat-'on  (26)  impossible. 
(See  figure  on  next  page). 

Hovfng  along  the  tangent  In  the  direction  of 
increasesing  T  (or  decreasing  V)  one  cones  to  lower  values 
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of  entropy,  whereas  we  already  know  that  in  the  case  Tg  >  T^ 
tl'G  point  2  should  correspond  to  a  higher  entropy  than  1.  But 
for   the  case  T.,  <  T,  the  more  refined  argument  using  the 
equation  (26)  seems  necessary. 


*  i  r  ',1  r>  5  T 


In  view  of  the  fact  that  T-,  -  D,    is  positive  or 
negative  accordnig  to  whether  T^  <  T,  ,  we  can  give  our  in- 
equality the  simple  form 

(D,  +  Dg  -  2T-^)-f'  (T^)  >  T^  -  D^  . 
Interchanging^  the  notations  T,   and  T^   one  finds  also 
(2T2  -  Dj  -  Dgj-fMTg)  <D2  -  Tg  . 


§ 


4.  Concluding  remarks.    The  third  consequence 


which  one  should  ho  able  to  drav/  from  the  condition  vi/f  >  0 

is  the  fact  that,  given  T,   and  T    there  is  s  single  interval 

of  the  T^  -  axis  (contained  In  the  interval  T  <  Tg  <  Tg  ) 

where  Q  =  1  -  4q-,q^  becomes  negative.   But  I  have  not  yet 

completely  succeeded  in  settling  this  point. 

Instead,  I  propose  to  show  by  a  quite  trivial 

argument  that  B(T2)  is  of  trough-like  shape,  more  exactly: 

COl^IDENTTAL 


28. 


Given  T,  ,  the  positive  function 

B  =  D^  4-  Lg  -  T^  -  Tg 

of     Tg      is   increasing  for     Tg  >  T^    ,    decreasing  for  Tg  <  T^    , 

Let  us    cotnplutc    the  derivative      '    with  respect   to   Tg, 
B'    =   Dj   +   D^   -   1    . 

The   equations 

D^Dg   =   T^Tg    ,  ^2   -    ^1      ^   ^(^2^    -    ^(^1' 

yield 

^P2  "  ^2  ^^1  =  ^1  '    ^2-4  =g'(Tg), 
consequently 

D»  (D^  +  Lg)   =  T^   -  D^g'(Tg)  . 


^^/ 


'•^1   '-    \^      =    -^1   -*■  ^2g'(^2^ 


(Dj  +  D^  -  1)(D^^  +  Dg)   =   (2T^  -  D^  -  Dg)  + 

+  (Dg  -  D^)(l  +  2f'(T2))   =   2[(T^  -  L\  )  +  (Dg  ~  D^)f'(T2)] 

This  formula  puts  the  sign  of  the  derivative  in  evidence. 


CONFIDENTIAL 


Date 

Due 

an  13  u 

f 

.-.--^      y 

W  •',>,, 

\ 

■/^ 

f-vV^'  '"^  ■ 

\   iC)7lR 

^^^ 

i  (JiJN  - 

\ 

.  A. . .    «•    n    ^1 

^fl'' 

JftKlGi 

JO"' 

D 

PRINTED 

IN  U.  S.  A. 

NYU 

AMG- 

18 

Weyl 

A   schema   ■P'— 


c.l 


N.  Y.  U.  Institute  of 
Mathematical  Sciences 

25  Waverly  Place 
New  York  3,  N.  Y. 


